Radiative decays of light vector mesons in a quark level linear sigma model 
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I. INTRODUCTION 



M. Napsuciale, S. Rodriguez and E. Alvarado-Anell 
Instituto de Fisica, Universidad de Guanajuato, 
Lomas del Bosque 103, Pracc. Lomas del Campestre 
37150, Leon, Guanajuato, Mexico 

We calculate the — > 77, P^"f and V'^ — > V''^77 decays in the framework of a C/(3) ® 

f/(3) linear sigma model which includes constituent quarks. For the first two decays this approach 
improves results based on the anomalous Wess-Zumino term, with contributions due to SU{3) 
symmetry breaking and vector mixing. The (j) — > {ui, p)^'f decays are dominated by resonant rj' 
exchange . Our calculation for the later decays improves and update similar calculations in the 
-closely related- framework of vector meson dominance. We obtain BR{cj> pj"f) ~ 2.5 x fO^^ and 
BR{(I} ^ ujjj) = 2.8 X 10"*^ within the scope of the high-luminosity (p factories. 
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10 ■ Radiative decays of mesons have been a useful tool in the past in the underlying of the dynamics of non-perturbative 
QCD. In spite of being one of the earhest tests for the internal structure of hadrons, the interest on this kind of 
processes is by no means exhausted. Indeed, there exists renewed interest on vector meson radiative decays [1] due 
to the possibihties opened by the high luminosity e~^e~ machines which are allowing us to study rare decay modes of 
QQ vector mesons. In particular many modes have been recently measured by the experimental groups at Novosibirsk [2] 
and KLOE Coll. at Frascati (j> factory started to release even more accurate data [3,4]. Hopefully this data will help 
us in discriminating the different proposals for the dynamics arising from QCD which governs the different channels. 

In the following we study the specific processes P° 77, V'~' P*'7 and V'^'yy decays. In QCD the 

P'^ 77 and P'^j decays are induced by the contributions of external vector currents to the axial anomalies 

which can be summarized in the Wess-Zumino term [5]. The V'^ y*^77 decays were previously calculated [6] in 
the framework of heavy vector chiral perturbation theory (HVCHPT) [7]. In this theory the V'^ V'^jj decays 
I '' are induced by loops of charged pseudoscalar mesons and the corresponding branching ratios turn out to be ratherly 
^H. small (w 10^9). The importance of intermediate pseudoscalar contributions to the same processes was noticed in 
(— I ' [8]. In this picture the V'^ — > V'^jj decays go through the chain V'^ —>■ P°7 V'^-j-j where the anomalous decays 
• V'^ — !■ P°7 play a prominent rol. The most important effect here is the possibility of resonant contributions which 
. !^ I enhance the corresponding branching ratios and in the case of — > P77 brings it close to the upper limit encountered 
' by the CMD-2 Collaboration [9]. Indeed, the branching ratio calculated in [8] for this decay is ,B = 1.3 x 10"'' to be 
Jh , compared with the upper bound reported in [9] Bexp < 5 x 10"''. The possibility of improving this bound at DA$NE 
■ - - 1 makes worthy to reconsider these decays. 

In this work we study the P° 77, V'^ P'^7 and V'^ — > y°77 decays in the framework of a U{3) (X) U{3) linear 
sigma model comprising Ua{^) symmetry breaking due to strong dynamics and including constituent quarks. This is 
a generalization to the U{3) (g) U{3) version of the U{2) (g) U{2) model used in [10] to calculate the cj — > /jtt transition 
and the uj —^ Sir decay. In this model the contributions to the axial anomalies due to external vector fields can be 
explicitly calculated yielding sharp predictions for the P° 77 and y° P°7 decays since we work with the physical 
meson fields. The advantage of the explicit introduction of fermion fields is that, on the one hand it makes possible 
to study effects due to SU{3) symmetry breaking and on the other hand it also allow us to study the contribution of 
resonant intermediate pseudoscalar states to more involved decays such as —> V'^jj. 

II. (7(3) x C/(3) LINEAR SIGMA MODEL WITH QUARKS (QLaM) 

The generalization of the model presented in [10] to its U (3)(XiC/(3) version can be done by studying the C7L(3)(Xif/fl(3) 
invariants constructed from B ~ S + iP. These invariants were studied long ago [11] 

X = tr{BB'i), Y^tr{BB'^f, Z ^ detB + detB^ . (1) 
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The first two terms arc ?7L(3)(Xit/i^(3) invariant whereas the Z term explicitly breaks the symmetry but respects 

the SUa^) X t^y(3) symmetry. The generalization of the Lagrangian presented in [10] is 



L = q[iD^r + V2g{S + ij5P)]q + ^triD^BD^'B^) - ^tr{BB^) 
^tr{BB^f - j{tr{BB^)f + tr{CS) - p{detB + detB^). 



(2) 



where the scalar and pseudoescalar fields are given by 5 = -^A*s% P = -^^^P^, with A* i = 1...7 the conventional 

Gell-Mann matrices and we use the flavor matrices Ans = -Dio£((l, 1, 0) and Ag = V2Diag (0,0,1). Explicitly, the 
scalar, pseudoscalar and vector fields are 
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Vector mesons are introduced as external fields through the covariant derivatives 



D^q= (d^ + i-^V^]q D^B = d^B + t^[V^,B]. 
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(4) 



The mesonic sector of this theory has been analyzed in detail in [12-14]. In particular, the vacuum expectation values 
(vevs) of scalars a = ^^'^ , b=< as > are related to the pseudoscalar decay constants as 



U = V2a fk=-j={a + b). 



(5) 



As for the fermionic sector, quarks acquire mass due to the spontaneous breaking of chiral symmetry. The constituent 
quark mass matrix generated by this mechanism is Mg = y/2g < S >, where < S >= V = diag{a, a, b). Explicitly 



mu=md = V^ga = gf,, 
rus = \pigb = g{2fk - U). 



(6) 



In the calculations of the physical processes to be considered below, the free parameters g, iriq, mg (q = u, d) enter in 
the combinations g/rUq, g/rus which are determined by the weak decay constants of pseudoscalars as 
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(7) 



The free parameter g^ appearing in the covariant derivative (4) can be fixed invoking vector meson dominance of the 
electromagnetic form factors of constituent quarks, i.e. VMD implemented at the constituent quark level (QVMD), 
or directly from the p tt+tt" decay. QVMD hypothesis was shown to work fairly well in the U{2) ® U{2) version 
of the model [10]. Under this hypothesis, the isovector contribution to the qq^ interaction is mediated by the rho 
meson. The corresponding amplitude is 



JvlVMD — e /p7 - 



2 2 I V"7i/W 

q-' — mj, 2 



(8) 
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where e stands for the photon polarization vector and fp^ denotes the p — 7 couphng. On the other hand, the amph- 
tude for the ggj transition can be characterized on the base of Lorentz covariance, gauge invariance, parity etc. as 



M 



qu(T'' 



(9) 



where stands for the quark charge and Fi{q'^), F2{q^) are form factors which can not be fixed on symmetry 
arguments alone. Eq. (8) has been written such that Fi{q^) is normalized to -F'i(O) = 1. In the soft limit {q^ 0) 
only the charge form factor contributes and comparing with Eq.(8) we obtain 



fp 



(10) 



It still remain to fix fpj. This parameter can be extracted from the measured p ^ e+e" decay. Using the central value 
for the measured branching ratio for this channel [18] we obtain Qy = 5.0. As a consistency check for this estimate 
wc can evaluate the isoscalar contribution to the qq'y interaction mediated by the ui meson. A similar analysis for the 
isoscalar form factor yields 



gv 



(11) 



Extracting the 0^ — 7 coupling from the u) — > e+e" decay we obtain = 5.6, close to the previous estimate. A 
calculation for the electromagnetic form factor of the s quark yields a similar value (gf„ = 6.0). Finally the coupling 
gy can also be directly extracted from the p° tt+tt" decay. From the Lagrangian (2) we get gp.^^ = g^. Using 
the measured branching ratio for tt+tt" we obtain g^ = 6.0, fairly consistent with the former estimates. In the 

whole, the vector coupling lies in the range € [5.0,6.0]. We will use gy = 5.5 ± 0.5 in the numerical calculations 
below. 



III. P ^ 77 AND V P-y DECAYS 

A. P°(p) 7(^1, £1)7(^2, £2) decays. 

The P° 77 decays are induced by quark loops in this model. These contributions are finite and we obtain sharp 
predictions for these decays. The invariant amplitude for the — > 77 decay as calculated in the model is 



M{Tr° 77) = C^77e(ei,e2,fci,/C2) /(p^) 
where e{p, q, r, I) is a shorthand notation for Cfj^vpaP^q'^r^l'^ and 

C^77 = 8imqgNc{el - el), 



(12) 
(13) 



with Nc the number of colors. The loop integral is given by 



1 /-l-X 



dy 



1 — x(l — x)^ + xy^ — ie 



(14) 



In the chiral limit {p^ = = 0) this integral yields 



/ = 



327r2m2 ' 



(15) 



and the invariant amplitude in the chiral limit reduces to 
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A^(7r° ^ 77) = -^e(ei,£2,fci,fc2) (16) 

where we used — = 4-- This resuh coincides with the amplitude calculated from the Wess-Zumino term ^. This 
amplitude yields the decay width 

r(7rO ^ ll)QLaM = = 7.63 eV (17) 

to be compared with the experimental result r(7r'^ — > 77)exp = (7.79 ± 0.56) eV. 

As for the 7] ^ 77 and rj' 77 wc must take care of the mixing of the 77ns and r/g in this case. The amplitudes for 
the physical processes are related to those for the flavor fields as follows 

Mri-y-y = A^))„377 COS (j)p — A4ri^^i~, sin d>p , , 

Mri'-y^ — Mri„sn ''i^ 't'p + -^Vsll COS (pp. 

The calculation of the amplitude for -q^s 77 goes along the lines of the case of the 7r°. We obtain 

■Mrinsji = (^r,ns'ri^{^u^2,ki,k2)I{ml^J (19) 

where 

Cn„.j-r = &imqgN,{el + el). (20) 
In the chiral limit / is given by (15) and the invariant amplitude reads 

5ct 

Miiins ^ -yy) = -—^e{si,S2,ki,k2) (21) 
where we used Eq. (7). Similar calculations for the r/s 77 transition yield 

■A4(?7s ^ 77) = ., f/^" , . e{ei,S2,ki,k2) (22) 
Finally, the invariant amplitudes for the physical processes ^ 77 y ??' ^ 77 are 

■^')77 = (Scos^p - sin (j)p^ e(ei, £2,^1,^2) ^^^^ 

Mn'j-f = 3^ (bsincpp + 2/f-/„ ^p) e(£i, £2, fci, ^2) 

Notice that these amplitudes reproduce the results of an effective theory of mesons in which the contributions of the 

electromagnetic external sources to the axial anomalies are introduced using 't Hooft's anomaly matching conditions 
[15]. As discussed in [15] from Eqs.(23) we obtain 0p S [38.4°, 41.0°], consistent with other estimates [16] and with 
recent measurements by the KLOE Coll. ^^^f = 41.8tJ J [3]. 

B. V{q,r]) — » P{p)'y{k,e) decays 

The most general form of the invariant amplitude for the V{q,r]) P{p)^{k,e) decays is dictated by Lorentz 
covariance, gauge invariance, parity etc. as 

M = Aype{e, rj, k, q) (24) 
where A^p is a form factor. The decay widths, written in terms of these form factors are 



^This is hardly a surprising result since the anomaly depend on the fermion content of the theory. The anomaly matching 
condition by 't Hooft used e.g. in [15], is satisfied here due to the same fermionic content in the model and in QCD. 



4 



T{V ^ P7) = 



967r 



1 - 



3 



(25) 



The value of these form factors on the mass-shell can be extracted from the existing measurements for the cor- 
responding decay widths against which the predictions of particular models must be compared. In the QLaM the 
P^j decays are induced by quark loops. The importance of constituent quark loops in the description of PVV 
vertices was stressed in [17]. The calculation of the amplitude for the p° — > w^j decay in the QLaM yields 



(26) 



where 



= 4m,<7.5iVe(e„ + ed)/(m2). 



(27) 



In this case the loop integral contains spurious imaginary parts coming from unphysical on-shell qq channel in the 
loops. This is due to the lack of a confining mechanism for the constituent quarks. We approximate this integral to 
its value at = and assume that this value is not too different from its value at the physical point in the presence 
of confinement. This is a reasonable assumption which, as we shall see below, is supported by experimental data. 
Under this assumption the value predicted by the QLcrM for Ap^ is 



aQLctM _ ^9v 



(28) 



In a similar way we calculate the kinematicaly allowed V P°7 and V^j decays, taking care of the mixing 

between the Vas and to form the physical <!) and u). Our results are shown in Table I where we defined 
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(29) 



We used gy = 5.5±0.5 as estimated previously, U = 91.9±3.5 MeV, fk = 112.9±1.27Mey [18] and cfy^ = (39.7±1.3)° 
in the numerical calculations. The vector mixing angle 0,,, is extracted from the experimental value for 7r*^7 
which yields (j)y = (3.31 ±0.42)°. The theoretical uncertainties come mainly from these quantities. The experimental 
amplitudes are obtained from the corresponding branching ratios [18]. 



Decay 


Amplitude(j4„p) 




A?/"^*^(GeT/-i) 


\A%^^\{GeV-') 




A 




0.229 ± 0.022 


0.294 ± 0.037 


p° Vl 


Bcos4>p 




0.529 ± 0.053 


0.562 ± 0.051 


u! — » 7r°7 


Bcos4>v 




0.687 ± 0.067 


0.711 ±0.016 


a; — » 777 


Acos(j)pCos4>v + Csin(f)pSinq 




0.164 ±0.017 


0.161 ±0.013 


V' - P°l 


Bsin(j)p 




0.439 ± 0.045 


0.388 ±0.006 


f] 0)7 


Asin4>pCos4>v — Ccos4>pSinci 




0.160 ±0.015 


0.137 ±0.007 


(p 7r°7* 


Bsirupv 




0.039 ± 0.001 


0.039 ± 0.001 


^ r?7 


Asin4>vCos4>p — Ccos(j)„sinci 


'>P 


0.211 ± 0.020 


0.209 ±0.002 




Asin(j)vSin4>p + Ccos(f)vCosci 




-0.233 ± 0.023 


0.225 ± 0.025 



TABLE I. Amplitudes for the V — ► P7 and P — ► V'l decays. The process marked as (*) is used as input to fix the vector mixing angle. 
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IV. V DECAYS 

The V — > ^'77 decays were firstly studied in [6] in the context of Heavy Vector Chiral Perturbation Theory 
(HVCHPT) [7]. There are two mechanisms contributing to those decays in this formahsm. The first one is the decay 
chain V VP ^'77. The second possibility is through loops of charged mesons. Furthermore, it was shown in 
[6] that these mechanisms do not interfere in the spin averaged decay rates. 

On the experimental side, the decay cfi p^^ has been already tested by the CMD-2 Coll. as a byproduct of their 
analysis of the ^ 7r"'"7r~7r° decay [9]. Although they are not able to discriminate if the later decay proceeds through 
the chain </> — » p-n^ — > 7r+7r~7r'^ or it is a direct decay, they can distinguish events where the two final photons come 
from the -k^ and pose an upper bound to the branching ratio for (j) — > P77 proceeding through other mechanisms. 
They obtain 

BR{(I) ^ P77)exp < 5 X 10-^ (30) 

In the context of HVCHPT, in addition to the (j) P'^^ Pll contributions which have been removed in Eq. [30] , 
this decay can proceed only through loops of charged mesons. The corresponding branching ratios were calculated in 
[6] as 



BR{4> ^ p^^)hvchpt = 5.8 x 10-^( 0%)^ . 
BR{ct> -> u:^^)hvchpt = 4.2 x > 



where the low energy constant is related to the VV'P couplings. Estimates for this constant using constituent 
quark model and chiral quark model yield 92 ~ 1 thus predicting extremely small branching ratios. 

The {p, a;)77 decays were reconsidered in [8] in the framework of Vector Meson Dominance. In this frame- 
work, in addition to the possibilities appearing in HVCHPT, the V — > ^'77 decay can proceed through the chain 
V P-f — > y77. The contributions from intermediate pseudoscalar states enhance the corresponding branching 
ratios to [8] 



BR{(I) p-fj)vMD « 1.3 X 10- 
BR{(I) ojjj)vMD ~ 1.5 X 10" 



.5 , (32) 



several orders of magnitude larger than those predicted by HVCHPT. The contributions of intermediate pseudoscalars 
do not appear in HVCHPT due to the conservation of heavy vector meson number inherent to this theory. 

Let us consider these decays in the framework of QLcrM. In this formalism the V V'"f"f are induced by quark loops. 
As discussed in the previous section, the assumed low variation of the form factors from — ^ rriy is supported by 
the experimental data onV P°7 and P° — » Vj decays. Here we will work with the same assumptions. Under these 
circumstances, our calculation is very similar to that performed in [8]. However, there are some differences which 
turn out to be relevant in the numerics. Firstly, our formalism incorporates SU (3) symmetry breaking, secondly we 
take into account mixing of vector mesons. The (/> {oj,p)jj decays are dominated by the exchange of on-shell 77'. 
Although there are also contributions coming from tt and rj exchange, the (j) 7r°7 decay is suppressed by the OZI 
rule (see Table I) and rj contributions turn out to be negligible as compared with the on-shell r]' exchange and will 
not be considered here. 

The invariant amplitude for <j){qi,r]i) p(92) ??2)7(^i) £1)7(^2, £2) is 



where 



with 



M = /te(?7i,£i,fci,gi)e(?72,£2,A;2,92) + /■ue(??i, £2, ^2, 9i)e(?72, £1, A;i, 92) (33) 



A 1 , A , A 1 , A , 

/* — ~7 ^2~ /w — ^2~ \^^) 



^ - iVn'rUr,', s = (91 - ^2)^, t={qi-kif, u={qi-k2f. (35) 
This invariant amplitude corresponds to t]' exchange in t and u channels. The average squared amplitude is 
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\Mf = ^ [\ftf{t - M^fit - m^f + \Uf{u - M^fiu - n?f + Re{ftf:)F{s,t, u)] (36) 

where 

F{s,t,u) = s'^M^m^ + {M^m^ -tuf (37) 

which is cxphcitly invariant under t ^ u as required by Bose symmetry. Here, M, m stand for the masses of (j) and 
the final vector meson respectively. The decay width is given by 



dr 1 11'"* 



ds (27r)3 32M3 2 

where 



\Mfdt (38) 

to 



to,i = |[(M2 + m^-s)T Vi^^ + m2 - s)2 - Am'^M^]. (39) 

If we integrate numerically this equation using the experimental values for the (f> rj'j, r)' p^f and r/ uj"f 
couplings listed in Table I, and the values for the masses quoted in [18] we obtain the following branching ratios 

BR{(I) P77) = 2.3 X 10-^ BR{(t) W77) = 2.6 x 10"^ (40) 

If instead we use the values for the couplings as predicted by the model also listed in Table I we get 

BR{(I, ^ p77)QiaM = 2.5 X 10-^ BR{(f) ^ u}'y'y)QLaM = 2.8 x 10"^ (41) 

These results are roughly five times smaller than those reported in [8]. It is instructive to understand where the 
differences come from. To this end let us analyze the (prj'^ coupling in both schemes. In our formalism, this coupling 
is given by 

9<t>ri'i = 0^2 J sin(j)ysin(j)p - — —cos(j)yCOS(j)p, (42) 



whereas the coupling used in [8] , when we use (pp = + a where sina = y/2/3, cosa = •\/l/3, can be written as 

94>ri'j = 4^2^ cx>S(l)p. (43) 

Clearly, modulo an irrelevant global sign, our analytical results reduce to those presented in [8] when we take the 
SU(3) limit (fk = /tt) and consider ideal mixing for vector mesons 0„ = 0. However, even if we take (f>y = 0, the SU(3) 
corrections introduce a factor ( ^^''^^^ )^ = (1.44)^ — 2.04 in the branching ratios. Further sources of discrepancy due 
to SU(3)symmetry breaking enter in the value used in the numerics for the pseudoscalar mixing angle. The extraction 
of this angle from 77 77 and 77' 77 is also sensitive to SU(3) symmetry breaking as can be seen in Eq.[23]. Our 
result, 4>p = 39.7 ± 1.3, is consistent with the recent measurement by the KLOE Coll. (1)'^^ = 41.8^16 ^^'^ larger 
than the one used in [8], namely cpp = 9 + 54.7° — 34.7°. Additional sources of discrepancy come from the non-ideal 
vector mixing angle and the numerical value used for g^,. The importance of all these effects is clearly exhibited in the 
calculation of the branching ratio for the ([) rj'^ decay. Indeed, this branching ratio was calculated to be 2.1 x 10^^ 
in [8]. In our formalism we obtain 6.8 x 10^^ to be compared with the recent measurement by the KLOE Coll. [3], 
namely BR{(j) 'ri'j)exp = (6.1 ± 0.61 ± 0.43) x 10~^. In summary, SU(3) corrections are necessary for the correct 
description of </>?7'7, hence they are also relevant for (j) (p, 0^)77 which are dominated by on-shell rj' exchange. For 
the latter decays additional sources for the different numerics come from the rj' and rj'ujj couplings, whose ratio 
deviates from the value of 3 predicted by SU(3)(see Table I), and from phase space integral. 

In any case, the (j) P77 decay is still close to the upper bound in Eq.(30) and both decays are within the reach 
of DA$NE. 



V. CONCLUSIONS 



In this work we study the decays P° 77, y° P°7 and V° V'^jj, in the framework of an ^(3) ® U (3) linear 
sigma model which includes constituent quarks. All these decays are induced by quark loops in this model. For the 
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P'^ —> 77 decays wc reproduce results previously obtained considering the contributions of external electromagnetic 
fields to the axial anomalies in the framework of the C/(3) (8) U{3) linear sigma model and using 't Hooft's anomaly 
matching conditions [15]. In the case of — > P^j decays there are spurious contributions coming from the opening 
of on-shell qq pairs in the loops due to the lack of a mechanism for confinement. Wc avoid this problem assuming 
that the corresponding form factors are smoothly varying functions of q^ and calculating their value at q^ = 0. 
Comparison with experimental results shown in Table I support this picture. The calculations for — > V'^'y-y under 
these approximations improve those performed in [8] by incorporating effects due to SU (3) symmetry breaking and 
the mixing of vector mesons. Our numerical results turn out to be smaller than those obtained in [8] roughly by 
a factor of 5. This difference can be traced to SU{3) symmetry breaking, non-ideal mixing for vector mesons and 
numerical input -which manifests in their large prediction for the branching ratio of the (/) rj'"/ decay- and also 
to phase space numerical integration. Our predictions for branching ratios of the most interesting processes are: 
BR{(j) — > p^^)QLaM = 2.5 X 10~^, BR{(j) — > u)j'y)QLcrM = 2.8 X 10~^. In the former case the calculated BR is quite 
close to the upper bound obtained by the CMD-2 Coll. BR{(f) pj^)exp < 5 x 10~^ [9] and these decays could be 
detected for the first time at DA$NE. 
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